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Abstract 

We propose a model for supernovae la explosions based on a phase transition to a supersymmetric 
state which becomes the active trigger for the deflagration starting the explosion in an isolated sub- 
Chandrasekhar white dwarf star. With two free parameters we fit the rate and several properties 
of type la supernovae and address the gap in the supermassive black hole mass distribution. One 
parameter is a critical density fit to about 3 • 10 7 g/cc while the other has the units of a space time 
volume and is found to be of order 0.05 Gyr R E where Re is the earth radius. The model involves 
a phase transition to an exact supersymmetry in a small core of a dense star. 
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I. INTRODUCTION 



It is central to string landscape ideas that the universe can exist in a (large) number of 
local minima in an effective potential and makes transitions between them. One of these is 
perhaps a local minimum with a very large vacuum energy describing the inflationary era. 
A second one is our broken supersymmetric (susy) world described perhaps as the Minimal 
Supersymmetric Standard Model (MSSM) with soft breaking parameters lifting the masses 
of the susy partners to hundreds of GeV and with a small positive vacuum energy. It is likely 
that the world of exact susy described by the MSSM with soft breaking parameters set to 
zero and vanishing vacuum energy also exists. In this world the quarks and leptons would 
have degenerate spin zero partners as would their baryon and meson bound states. Such 
an exactly supersymmetric universe could be the ground state of the string landscape or at 
least a long lived intermediate state with the true ground state being a state of negative 
vacuum energy as in the conformal field theory/ Anti-deSitter correspondence (CFT/AdS). 
The transition between our universe and the zero vacuum energy state is describable as a 
decay of the false vacuum. The transition probability would presumably be affected by the 
presence of matter and we assume that it is enhanced rather than suppressed as suggested 
by analyses in lower dimensional models and by the fact that the vacuum decay probability 
is an increasing function of the energy density difference between the two states. 

Supernova la (SN la) explosions have been successfully used to provide strong evidence 
that the expansion of the universe is accelerating. Vital to the argument is the understand- 
ing of the calibration of the light curves, and for that we need insight into the explosion 
mechanism. The primary standard model picture for SN la prior to 2010 was the single 
degenerate model [lj whereby a white dwarf accretes matter from a main sequence binary 
partner until it approaches the Chandrasekhar limit and collapses. This model has been 
greatly disfavored by recent x-ray analysis Although many doubly degenerate models 
whereby two white dwarfs merge or collide are perhaps still viable, such models involve a 
host of free parameters with, consequently, little predictive power. It is not known how many 
suitable binary systems exist or what would happen to the many systems without the fine- 
tuned orbital parameters needed to reproduce observations. More complicated, shrouded, 
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singly degenerate accretion models could also be devised to circumvent the result of [2J, but 
these add free parameters and extra complexity. 

In view of the doubts thrown into the debate of what SN la really are, the time may be 
ripe to reconsider phase transition models in parallel with further exploration of standard 



model possibilities. In a recent review sub-Chandrasekhar explosions, of which our 
current model is an example, are identified as one of two scenarios which might ultimately 
explain the majority of SN la events. 

The primary features of SN la which we address in our two parameter model based on 
string landscape and vacuum decay ideas are: 

1. no hydrogen in emission; 

2. homogeneity; 

3. light curve dominated by 56 Ni production and decay; 

4. rate: About one per century per galaxy 

(more in galaxies with enhanced star production but also observed in old galaxies); 

5. some evidence for appreciable departures from spherical symmetry; 

Item 1 is naturally fit if the SN la progenitors are white dwarf stars with no circumstellar 
hydrogen as occurs both in our model and in doubly degenerate accretion models. 

As discussed below, our model predicts a narrow range of progenitor white dwarf masses 
which naturally results in the SN la homogeneity (Item 2 above) with, however, some spread 
in energy release. 

The supersymmetric (susy) bubble grows in the present model from a sub-atomic size 
to a macroscopic size which, however, is limited by the lack of degeneracy pressure in the 
susy phase. The energy release induces fusion in the outside normal phase leading to the 
observed 56 Ni production which is a natural endpoint of fusion in a carbon-oxygen white 
dwarf (Item 3). 

The rate of SN la is easily fit in our model using the known mass spectrum of white 
dwarf stars (Item 4). The observed rate yields a strong relation between the model's two 
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free parameters. This feeds back into the narrowness prediction and also naturally yields 
the lifetime observations of item 4. 

Although the nucleation of the susy bubble begins most probably at the stellar center, the 
average location of the nucleation point is quantitatively predicted to be off-center leading 
naturally to an off-center explosion (item 5). This average location is predicted to be just 
inside the sphere containing the entire region above critical density with a broad distribution 
throughout this region. 

A fit to one of our two free parameters suggests that the supersymmetric (susy) phase 
transition described above occurs in single white dwarfs of mass from about 0.9 to about 1.3 
solar masses, starting a deflagration front, which then triggers the explosion as a supernova 
ll, |4j . Since the data clearly show that Supernovae la occur both in stellar systems, which 
are constituted out of mostly old stars (elliptical galaxies) , and in stellar systems made of old 
and many new stars (late Hubble type spirals) [5j, there has to be a very broad delay time 
distribution from making a white dwarf to its explosion. We emphasize that what we call 
here a delay time is more akin to the decay of a radioactive nucleus, so any particular white 
dwarf passing the threshold may still have a chance to go susy somewhat earlier and also 
somewhat later than the nominal delay time. For single stars a broad delay time distribution 
between creation and explosion constitutes a challenge, but a susy phase transition offers 
such a prediction. Super-massive black holes have been observed over a wide mass range, 
:'rom about 3 • 10 6 M Q to 3 • 1O 9 M , with some outliers to lower masses such as 3 • 10 5 M 
6] . There is also strong evidence for stellar black holes, in the approximate mass range 3-5 
Mq. However, there is currently no convincing evidence for a large number of black holes in 
the intermediate mass range 7[. These are either difficult to detect, or just may not exist. 
Observational arguments for the existence of a few do exist js]; expected are many. Here we 
pursue the notion, that such black holes can derive from the agglomeration of very massive 
stars: we explore the concept that across the mass range of the gap the final collapse of the 
star leads to a susy phase transition in the core. This leads to an additional energy input 
possibly pushing these stars into an enhanced mass loss during collapse, and may result in 
a final compact object of much less than 1O 5 M . 

From many points of view, the theory of violent astrophysical events would be less prob- 



1 



lematic if there were a new source of energy release beyond the standard model. For instance, 
in the standard model the mass of a white dwarf grows by accretion to the Chandrasekhar 
limit. However a star near this limit has a binding energy greater than available from fusion 
reactions making it impossible to totally unbind the star without borrowing energy from 
the gravitational energy released during the accretion. This energy must be repaid as the 
star explosively expands making the timing critically sensitive. Similarly in earlier years 
it was thought that the energy required to blast off the outer shell of massive stars in a 
Type II supernova was provided by neutrinos but many detailed studies have found that 
such supernovae would stall because of insufficient numbers of neutrinos and the weakness 



of neutrino interactions 



9|. In the case of the Type Ib/c supernovae of very massive stars 



magnetic fields may help cause the explosion 10l. Ill|. 



A phase transition to exact susy in dense stars could be nature's way to release the energy 
stored in a Pauli tower of fermions. A susy phase transition in white dwarf stars (WD's) 
could eliminate the need for accretion and allow an isolated white dwarf to explode with the 
required extra energy input and a correctly predicted rate. 

II. THE DECAY ASSUMPTION 

In the vacuum the probability per unit time per unit volume for the decay of the false 



vacuum is governed by the Coleman-DeLuccia formula 121 ] . 



P 

A c e~ B ^ (1) 



dtd 3 r 
with 

\ 21tx 2 S 4 

B{vac) = 2K^ (2) 

Aq" = To Vq is a parameter with the units of a space-time volume. S is the surface tension 
of the true vacuum bubble in the dominantly false vacuum background. In the case of a 
susy phase transition e is the difference between the vacuum energy density of our universe 
and the zero vacuum energy density of an exactly supersymmetric background. 



In dense matter it is expected that the transition is accelerated since the energy density 
difference between a broken susy ground state and that of an exact susy ground state with 
the same additive quantum numbers (baryon and lepton numbers) is dominated by the 
energy stored in the Pauli towers. Thus it is plausible that the phase transition formula 
could, to at least a first approximation, be given by replacing e with the total energy density 
difference e + Ape 2 . At present however, a rigorous proof of the acceleration exists only in 
lower dimensions 

Other attempts to model matter induced or catalyzed vacuum decay, although different 
from the present proposal, can be found in ref. [ijj]. 



In exact susy , 



mechanism 
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, the degeneracy of bosons and fermions plus the availability of a conversion 



16[ from a pair of fermions to a pair of partner scalars implies that Ape 2 is 



the Pauli excitation energy density of the fermions. The Fermi gas model predicts for this 
excitation energy in a heavy nucleus of N neutrons and Z protons 

AE 1 /,2JV\ K / a ,2Z, 



Ap = p 



^5/3 , (±^\B/3 



AM n c 2 2 V A ' A 
3(9vr) 2 / 3 hp 



0.02p . (3) 



40 M n cR 

Here M n is the nucleon mass, Rq is the nuclear constant 1.2 fm, and the approximate 
value is given for a nucleus with N = Z = A/2. This energy release is about three times 
that of standard hydrogen fusion and some twenty times greater than the carbon or oxygen 
fusion which could provide the energy release in a standard model SN la. 

Thus the parameter controlling the exponential factor in the transition rate would then 

be 



We would then have 



27tt 2 5 4 

~ 2hc(e + Ape 2 ) 3 ' ( ^ 



H 2 P 1 

■ e~ B . (5) 



dtd 3 r TqVq 

This equation should be taken as the basic assumption of this paper with the previous 
discussion of vacuum decay and possible connections with string landscape ideas providing 
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some motivation and inspiration. Although the vacuum decay theory, based as it is on the 
WKB approximation, is subject to corrections, the current paper attempts to explore the 
consequences of Eq.flS]) as written with tqVq and S treated as two free parameters. We do 
not provide a theory of susy breaking although, when one is eventually found, it should 
predict the values of our free parameters. We use the observation that the vacuum energy 
of our broken susy universe is positive while, in exact flat space susy, the vacuum energy 
vanishes. An auxiliary assumption is that the degenerate mass of fermions and bosons is 
equal to the fermion mass in the broken susy world. This would be as in the MSSM or one 
of its extensions with the soft susy breaking parameters set to zero. 

With our two parameters we fit the SN la rate and explain the energy release in SN la 
(approximately) as well as the narrowness of the progenitor mass distribution and the time 
delay distribution from the white dwarf birth i.e. the onset of degeneracy. Note that white 
dwarfs in our predicted mass range do not spend appreciable time in the red giant phase. 

In dense matter containing heavy nuclei (above helium in which all nucleons are in the 
lowest level), e is negligible compared to Ape 2 . 

We would therefore write, using the Fermi gas model and replacing S in terms of a critical 
density p c , 

= e ( ~> . (6) 

dtcPr T V y ' 

The white dwarf density, p, as a function of radius is given as in fig. [T] depending on the 
total mass of the star. The phase transition probability per unit time per unit volume, 
Eq. ([6]), increases rapidly with p until p becomes of order p c at which point it saturates. 
For denser media, the transition rate is proportional to the volume. Thus, if we choose p c 
as near the white dwarf density, the transition rate in neutron stars will be at least eight 
orders of magnitude lower than that in white dwarfs. If there was a degenerate quark gluon 
plasma in the very early universe, one could consider whether the universe might have turned 
supersymmetric at that stage thus ruling out our model but, in fact, the space time volume 
of this phase would have been too small compared to tqVq to effect the phase transition with 
an appreciable probability. We would estimate the probability as roughly 10~ 10 . 

For present purposes we will ignore the additional Pauli energy stored in electronic states. 
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III. GROWTH OF THE SUSY CORE 



In dense matter the critical radius, generalized from the vacuum decay result would be 

Rc = -^— 2 • (7) 

e + Ape 2 

Bubbles nucleated with less than this critical radius will be immediately quenched. This 
follows from the vacuum decay result if, in dense matter, the relevant quantity is the total 
energy density difference between the two states and not just the vacuum energy difference 
as in the basic assumption of Eq. A bubble of critical (or greater) radius will grow as 
long as its radius exceeds the density dependent critical radius of Eq. (j7j). Replacing S by 
a critical density as discussed above leads to a critical radius of sub-atomic size in a white 
dwarf but still much larger than a nuclear radius. In a medium of uniform density (such as 
the vacuum) the bubble will grow without limit but a susy bubble nucleated within a dense 
star will not grow beyond the boundary of the star since outside the star the critical radius 
is of galactic scale. However, there is also a hydrodynamic limit on the growth of the bubble 
in dense matter due to the absence of degeneracy pressure. The net effect is that the susy 
core never grows beyond a small fraction of the host star. 

It is important to work out the susy phase transition front growth: A phase transition 
with energy input can be seen as analogous to an HII region ionization front; the ionization 
injects energy and there is a surface discontinuity, usually a shock front. Writing for the 
speed of sound ahead of the front, in the normal white dwarf matter, Cwi, and for the speed 



of sound behind it, in the susy matter C susy , the advance speed Ud is ([171]) given by what 
is called a D-type front: 



Un - ^ • (8) 

^ ^ susy 

Since very likely Cyjd/C susj < 1, it is probable that the advance speed is below the speed 
of sound, and so the susy core grows only slowly. In lower density regions, the speed of 
sound is considerably less so the bubble growth is further slowed. Particles that are light in 
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both phases can pass freely through the phase boundary and induce normal fusion in the 
region outside the core. Normal fusion reactions outside of the susy core produce the Nickel 
that governs the SN la light curves. Since, according to calculations (see [3]), the Nickel 
seems to be produced at intermediate r, it is an important feature of our model that the 



susy core remains small compared to the radius of the star 16]. 

From Eq. ((6]) we see that each star of mass M has a characteristic lifetime r(M) defined 

by 



dt T V J r V t(M) 

For very dense stars V c is approximately the volume of the star at or above the critical 
density. To consider the transition in a dense helium star one would have to replace Eq. (j3J) 
by the appropriate formula for an electronic rather than nuclear degeneracy. For a given 
susy bubble surface tension, the transition probability per unit time per unit volume would 
then be much less than for a white dwarf star. In the absence of new star production this 
would lead to a number N(M, t) surviving a time t after achieving a constant density defined 
by 

N{M, t) = N(M, 0) e"' /r(M) (10) 

The complement of the implied probability defines the likelihood that a white dwarf of 
mass M will explode in time less than t. The current model can be ruled out if a significant 
sample of white dwarfs can be found with ages long compared to their r(M) defined in 
Eq. (}9"1) . At present only low mass white dwarfs are known to have existed over gigayear 
time scales. 

New white dwarfs will be produced at a declining rate as the raw material in the galaxy 
is depleted. We can, for simplicity, parametrize this in terms of a star formation time scale, 
Tf, and a normalization, A(M): 

dN(M,t) A(M)_ e _ t/ ^ (n) 



dt Tjj, 

If we combine Eq. ([9]) with this rate of production of new white dwarfs we would predict a 
time dependent number of white dwarfs given by 

N(M,t) = A(M) / (M) r (e"*M M ) - e -*/**) (12) 
V ; V ; (t(M) - r f ) v ; v ; 



Equation Q defines a natural lifetime r(M) for any broken susy object containing heavy 
nuclei, but for objects far below the white dwarf density, such as our current sun, this lifetime 
is, by our choice of the TqVq parameter, effectively infinite, i.e. many orders of magnitude 
longer than the current age of the universe. 

Just as it is statistically possible for an elementary particle to decay in a time less than 
its natural lifetime, a white dwarf in our model could also produce a supernova in a time 
less than r(M) although r(M) would be the average life of an ensemble of such stars. The 
probability for a phase transition to occur at some time > t after reaching its classical steady 
state would be e~ t l T ( M \ Note that for high mass white dwarfs (M > 0.9 M Q ) the time spent 
in the red giant and previous phases is small compared to a gigayear. 

It has been suggested [3] that the critical radius in vacuo should be the galactic radius 
(« 4.7 ■ 10 22 cm). This would lead to a surface tension S = 8.9 • 10 13 erg/cm 2 and then to 
p c = 3 • 10 5 g/cc. An independent estimate from applying the susy phase transition model 
to gamma ray bursts [ljj] leads to comparable values from p c « 10 6 g/cc to p c ~ 3 • 10 7 g/cc. 
Thus from two very different starting points a critical density near the white dwarf density 
is suggested. As we shall see, the edge of the high mass black hole distribution also suggests 
interesting new physics near this density. 

Without loss of generality one can choose Vo to be the maximum V c for all white dwarf 
masses at a given p c . Then the parameter tq represents, for given p c , the minimum quantum 
mechanical lifetime of white dwarfs against nucleation of a susy core. The lifetime of any 
other star is inversely proportional to its critical volume, V c . Ultimately, a star with such a 
susy core must either be totally disrupted by the susy energy release or collapse into a black 
hole due to the absence of degeneracy pressure. The time over which a star can survive with 
a susy core is dependent on the rate of susy energy release. 

In a single star model for SN la we obviously have no need to push a star's mass up with 
time via accretion from a binary partner. What we propose is commonly referred to as as 
a sub-Chandrasekhar model: the explosion of a white dwarf below the Chandrasekhar mass 
limit. In our model this is triggered by the phase transition to a susy phase inside a very 
small core of the white dwarf. Thus this model makes no distinction between C+O white 
dwarfs and Ne-0 white dwarfs, except in so far as they exist in different white dwarf mass 
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FIG. 1: The density of three white dwarfs as a function of radius relative to the Earth radius 
ve- Dashed lines give the range of previously discussed critical densities. Masses of the three 
white dwarf examples are given in Mq. These Chandrasekhar density distributions go into the 
calculation of V c in Eq. (|9j). 



ranges. 

The explosion rate can then be written as 

dN SNIa 



dt 



WD' 



G(Pc 



where 



G( Pc 



dM- 



1 dN WD V c 



(13) 



(14) 



iV W D dM V 

Here the white dwarf mass distribution - which can be modelled over the relevant mass 
range using the Salpeter mass function and the lifetime of stars on the main sequence - 
is taken to follow the high mass tail of the Sloan white dwarf sample of 4621 stars with 



T e ff > 12000-fT 20|, |2l|. The suggested white dwarf mass distribution in the high mass tail 
going into Eq. (1161) varies as M -4,78 although the result is not greatly sensitive to the exact 
behavior since the relevant mass range is small. 
iVwn in Eq. fTTBl is the total number o: 



at near 10 10 [22|. Statistics suggest 



white dwarfs in the Milky Way galaxy, estimated 
23] that SN la are about as common as SNe from 
higher masses. Assuming that all white dwarfs above some particular zero age main se- 
quence (ZAMS) mass and below 9 M & become SN la, and that all stars above ZAMS mass 
9 Mq become other kinds of SNe, suggests 24j that this critical mass is about 5.5 M Q , 
corresponding to a white dwarf mass of about 0.9 M & again, approximately. This is the 
mass at which the delay time is still less than the Hubble time, so that we still can get a 
few white dwarfs exploding at that mass within the age of a Galaxy. 
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FIG. 2: The relative delay time distribution, t/tq, as a function of mass for critical densities of 10 6 , 
10 7 , 10 8 , and 6.3 -10 s g/cc. For an approprate choice of p c , the narrowness of the SN la distribution 
is naturally explained. We superimpose the histogram of the observed mass distribution of white 
dwarfs above an observed temperature of 12,000 K 20J. For p c > 10 7 g/cc and to = 0.5 Gyr only 
high mass white dwarfs undergo the susy transition with a lifetime less than the current age of the 
universe. 



The SN la rate is found to be a function o 
in a galaxy and the rate of star production 



25|: 



both the total number of white dwarf stars 



dSNIa 



A- 



dM. 



dt "10 10 M Q 
with A = 0.04+™ 16 and B = 2.6 ± 1.1. 



B 



dt 



10 W M @ /Gyr 



(15) 



The observation that the supernova rate is greatly enhanced in galaxies with elevated star 
production rates suggests that there is a component of the supernova rate at small delay 
times. The observation that supernovae continue to be produced in older galaxies suggests 
that there is also a component with longer delay times. Standard model approaches to 
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FIG. 3: The function G(p c ) which gives via Eq. (|13p the supernova rate. p c is given in g/cc. 




(Gyr) 



FIG. 4: The relative lifetime distributions for the critical density of 3.55 • 10 7 g/cc, as a function 
of r for ro = 0.5 Gyr. The evident scaling behavior in Eq. (]16p shows that the graph is a function 
of r/ro only. The distribution is sharply peaked at r = To- 

SNIa strive to be compatible with this observation leading perhaps to a bi-modal delay time 
hypothesis but no firm prediction. 

In our model there is a natural spike in the lifetime distribution near some minimum ro 
plus a tail to much higher lifetimes. This distribution is given by 
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As can be seen from fig. [2] r(M) is approximately parabolic on a log-linear scale. This 
implies that near t = tq there is an integrable singularity: 

dN 1 



(17) 



dr T ^/ln t/tq 

Without relying on this approximate behavior the numerical solution of Eq. (I16p is given 
in fig.H] showing clearly the spike at low lifetimes. On the y axis in fig. H] is the predicted 
number of white dwarfs with given lifetime in a sample of the size of the Sloan white dwarf 
sample of 4621 stars 20|, |21[. The logarithmic spike predicted in the present model is similar 
to and probably observationally indistinguishable from the empirical fit 

dN 1 

* " Tr (18) 
suggested in ref. 5[. These estimates suggest a critical density of order 3- 10 7 g/cc, which we 
have used in fig. HI The instantaneous SN la rate as given in Eq. (1131) is a function only of the 
total number of white dwarfs in the sensitive mass range but, during periods of elevated star 
production rate, this number is also temporarily elevated leading to the B term in Eq. f lT5|) . 
As the star production rate declines a dip should appear in the white dwarf mass distribution 
at the mass of minimum lifetime although the statistics are too meager at present to see this 
in our galaxy. Another prediction of the current model with our fitted value of p c is that 
the total mass ejected in the supernova should range from 1.0 M to the Chandrasekhar 
mass. In a binary accretion model the total mass ejected should be very close to or above 
the Chandrasekhar mass. Current measurements are not sufficiently precise to rule out or 
confirm either model. 

We can check how much energy is released by the phase transition of a very small part of 
the white dwarf: If we require that the susy energy release plus that from fusion outside the 
susy core is equal to the binding energy of the outside shell plus some 10 51 ergs of kinetic 
energy we obtain that the susy core is very much smaller than the stellar radius and the 
proportional mass in Nickel can easily vary by over an order of magnitude. It is also possible 
that the variability in the distance of the susy nucleation point from the center as shown 
below in fig. [5] contributes to the variability in the amount of Nickel produced. 

For SN la cores we need less than 10 -3,5 M ( 4j) for an initiating phase transition. Since 
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the susy phase transition provides more energy than even hydrogen burning, even less matter 
than in normal models is required to initiate the deflagration. 

One proof of our approach could be if these very small compact susy objects could be 
found remaining after the explosion. In some respects these objects might resemble black 
holes of anomalously small mass. 

Folding the white dwarf mass function with the delay time distribution allows on the one 
hand to obtain the SN la progenitor mass distribution (see fig.EJ), and also to obtain the SN 
la distribution as a function of delay time (fig.HJ). Most SN la are produced after a lifetime 
near To although there is a tail extending to higher lifetimes (see fig-H]). With a small To we 
naturally understand the greater rate of SN la in galaxies with an elevated star production 
rate and the lower rate in our galaxy. 

As greater statistics of high mass white dwarfs are accumulated, we would predict a dip 
in the white dwarf mass distribution at the minimum of the t/tq curve for chosen critical 
density in fig.EJ Since our approach yields a very narrow mass distribution for those white 
dwarfs that explode as SN la, the energy distribution is also very tight and so a use as a 
standard candle is justified. The apparent luminosity is, therefore, a good measure of its 
distance. 

Data also show 26| a relatively large number of SN la with a delay time of order 5-10 8 yrs. 



Combining these two points yields the critical density and the critical time scale. Therefore 
a clear prediction of this model is a relatively sharp sub-Chandrasekhar mass distribution 
(see, e.g., ref. j^)' peaking at the minimum of the delay time distribution. Thus our 
estimates for the key parameters are p c ~ 3 • 10 7 g/cc, To — 5 • 10 8 yrs, resulting in a typical 
progenitor mass near 1.2 Mq, mostly an original C-0 composition, and an average SN la 
rate per galaxy of about 1/100 yrs, see Eq. (1131) and the accompanying fig. El Considering 



the range of uncertainty in the observations 



23 



26(|), the time scale could be yet shorter, 



implying a smaller value of G(p c ) and so a higher critical density; similarly, the fraction of 
all SNe turning into SN la could also be smaller, pointing in the same direction. This might 
be verifiable with a few well observed SN la. 

Although the most probable nucleation point of the susy core is at the stellar center, from 
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M/M Q 

FIG. 5: The mean nucleation point of the susy bubble relative to the stellar center with rms 
deviation indicated based on a critical density of 10 7 g/cc. For high mass stars the rms upper limit 
is approximately the radius at which the density is equal to the critical density. The upper line 
gives the stellar radius as a function of white dwarf mass. Re is the Earth radius. 

Eq. we can predict the mean distance from the stellar center of the nucleation (see fig. [5]). 
Thus the model gives a clear quantitative prediction of the amount of "off-centeredness" of 
the stellar explosion. The fusion to 56 Ni is presumably centered in the most dense part of 
the star outside the susy core regardless of where the susy energy release originates. 



A. Massive stars and the Black Hole Gap 



Applying the susy transition idea to massive stars as they agglomerate 28J within a dark 
matter clump in the very early universe also leads to a phase transition that may be fast 
enough to slow the collapse and eject much of the mass of the star before it ultimately 
collapses as a black hole (BH) at much lower mass. This would then help to explain the 
BH mass distribution gap between about 30 and about a million solar masses, where we 
observe very few BHs. Going through the numbers 29( suggests 10 5 solar masses as the 
upper limit for this mechanism since this is the dividing line to achieve white dwarf density 
before collapsing into a BH at higher masses. We note that the general relativistic instability 



of very massive stars found by Appenzeller and Fricke is at slightly higher mass ! 
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The maximum average density that a stellar conglomeration of mass M can achieve before 
becoming a black hole is 



M 



/W = 1.83 • 10 6 [jf^j g/cc ■ 
The peak density is expected to be 54 times greater (Chandrasekhar 1939) 



Ppeak = 9.9 ■ io 7 (rfirgA s/ cc • ( 20 ) 

The numerical factor in Eq. ( Tl9|) is a typical white dwarf density pwd = J~]§r- Thus with 
a critical density in this range, the susy phase transition predicts an edge in the supermassive 
black hole distribution near 1O 5 M . Above this mass a conglomeration can become a black 
hole before the critical density for a phase transition is reached. Below this mass the susy 



transition will take over and disrupt the collapse into a black ho 
core which will eventually become a susy black hole is estimated 



e. The volume of the susy 



161 ] to be less than 0.1% of 



the original stellar volume. Thus the present model could help explain the gap in the black 



hole mass distribution 



BE 



This then leaves the range of about 10 5 to about 10 6 solar masses to make the first 
generation of super-massive BHs; these BHs can merge and then may help to explain the 
observed supermassive BH mass distribution 



B. Supernova Characteristics 

Returning to supernova properties, a clear prediction of our model is then the absolute 
rate of SN la, see Eq. ( I13p and Eq. (I14p in conjunction with fig.[3]and the relative statistics 
in terms of mass, energetics, chemical composition, of white dwarfs which explode and with 
what delay. A key prediction of our model is the shape of the delay time distribution, 
see fig.Hl Detailed observations and modelling should confirm that the total mass extends 
significantly below the Chandrasekhar limit, and that the original white dwarf was composed 
mostly of carbon and oxygen, since these are the stars in the preferred mass range. A further 
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prediction is the existence of very small mass compact susy remnants, observable in their 
passage through the interstellar medium due to radiation from accreting material. 
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